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Abstract 

We present a preliminary study of a new class of two-input cellular automata 
called eventually number- conserving cellular automata characterized by the 
property of evolving after a finite number of time steps to states whose num- 
ber of active sites remains constant. Eventually number-conserving cellular 
automata are models of open systems of interacting particles, that is, sys- 
tem of particles interacting with the external world, The particle aspect of 
eventually number- conserving cellular automata can be emphasized by the 
motion representation of the cellular automaton evolution rule. This new 
class of cellular automata contains, as strict subclasses, number-conserving 
cellular automata, monotone cellular automata, and cellular automata emu- 
lating number-conserving ones. Our main objective is to show that they are 
not what one might naively think they are. 

1 Introduction 

Whereas cellular automata have been widely used to model complex sys- 
tems in which the local character of the evolution rule plays an essential 
role [U 12 El HI E], the theory of cellular automata is still in its infancy and 
very few exact results — as, for example, the necessary and sufficient condi- 
tion for a cellular automaton rule, either deterministic or probabilistic, to be 
number-conserving — have been obtained [UJ [7| |H1 El UH UJ] Recently most 
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of the theoretical research carried out into this field is devoted to definition 
and characterization of new families of rules. In order to extend the fam- 
ily of number-conserving rules a number of papers on monotone rules (the 
set of number-conserving rules is a subset of monotone rules) have been pub- 
lished ^2] ■ The complexity of problems concerning characterizations of these 
properties may be illustrated by a recent result ^3] showing that monotony 
is decidable in dimension one but not in higher dimensions. This paper is 
devoted to the study of a new class of cellular automata called eventually 
number- conserving cellular automata defined by the property of evolving af- 
ter a finite number of time steps to states whose number of active sites re- 
mains constant. The characterization of this new class of cellular automata 
is challenging since both number-conserving and monotone rules are special 
cases of eventually number-conserving rules, and it is not, therefore, obvious 
that eventual number conservation is a decidable property. 

A one-dimensional cellular automaton (CA) is a discrete dynamical sys- 
tem, which may be defined as follows. Let s : Z x N h Q be a function 
satisfying the equation 

s(i,t + l) = f(s(i - r e ,t),s(i -r e + l,t), ...,s(i + r r ,t)), (1) 

for all i G Z and all t G N where Z denotes the set of all integers, N the 
set of nonnegative integers, and Q a finite set of states, usually equal to 
{0,1,2, ... ,q — 1}. s(i, t) represents the state of site i at time t, and the 
mapping / : Q r *+ r r+ 1 _ > Q i s the local CA evolution rule. The positive 
integers 77 and r r are, respectively, the left and right radii of the rule. In 
what follows, / will be referred to as an n-input rule, where n is the number 
rg + r r + 1 of arguments of /. Any map c : Z^ i— > Q is called a configuration. 
Following Wolfram to each rule / we assign a code number N(f) such 
that 

N(f) = E x ^ ■ ■ ■ > x n )g 9 "" lxi+9 "" 2x2+ - +9 ° x ". 

Two-state, three-input cellular automata (CAs) are referred to as elementary 
CAs. 

In this paper, we will only consider finite (or periodic) CAs, and replace 
the set Z by the set Z^, of integers modulo L. Any element of the set Q L will 
be called a cyclic configuration of length L. 
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A one- dimensional q-state n-input CA rule f is number-conserving if, 
for all cyclic configurations of length L > n, it satisfies 

, X2, ■ ■ ■ , x n) + f( x 2, x 3, ■ ■ ■ j x ni x n+l) + ' ' ' 

+ f( x L, Xi..., X„_ 2 , X„_l) = Xi + X 2 H h Xi. (2) 

It can be shown [7] that: 

one- dimensional q-state n-input CA rule f is number-conserving if, 
and only if, for all (x±, x 2 , ■ ■ ■ , x n ) G Q n , it satisfies 



n-l 



f(Xl,X 2 , 




k 



The purpose of this paper is to study eventually number-conserving CA 
rules. This new class of CA rules, which is an extension of the class of 
number-conserving CA rules, is defined as follows: 

A one- dimensional q-state n-input CA rule f is eventually number-con- 
serving (ENC) if, after a finite number of iterations of rule f , it satisfies 
condition (HJ) for all cyclic configurations of length L > n. 

While number-conserving CAs may be viewed as models of isolated sys- 
tems of interacting particles in which processes of annihilation or creation of 
particles are forbidden, ENC CAs are models of open systems of interacting 
particles exchanging particles with the external world such that at after a 
finite transient time the numbers of particles moving from and to the sys- 
tem exactly counterbalance each other. Numerical simulations show that the 
constant number of active sites in configurations of the limit set of an ENC 
CA depend not only upon the initial number of active sites but also on the 
detailed structure of the initial configuration. 

In the two following sections we consider classes of CAs that one might 
think coincide with ENC CAs. In the last section the motion representation 
of CA rules is introduced. More explicitly than the usual rule table, it shows 
how to represent particle motion either inside the system or between the 
system and the external world, stressing the fact that ENC CAs can be 
viewed as systems of interacting particles. 
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2 Monotone CA rules 



A one- dimensional q-state n-input CA rule f is number-nondecreasing if, 
for all cyclic configurations of length L > n, it satisfies 

, X2, ■ ■ ■ , X n _i, X n ) + f(x2, 

+ f(x L , X! . . . , X n - 2 , Xn-l) > Xi + X 2 ~\ h X L . (4) 

The g-state n-input CA rule C(f), defined by 

C(f)(x 1 ,x 2 , . ..,x n ) = f{q- 1 - Xt,q- 1 — x 2 ,...,q- 1 - x n ), 

is called the conjugate of rule /. This definition implies that the conjugate 
of a number-nondecreasing CA rule / is number-nonincreasing, that is, for 
all cyclic configurations of length L > n, it satisfies 

f{xi , X21 ■ ■ ■ , X n _i, X n ) + f[X2, X3, . . . , X n , X n ^_ij -\- 

+ f(x L , Xi..., X n -2, X n -i) < X! + X 2 H h X L . (5) 

A CA rule is monotone if it is either number-nondecreasing or number- 
nonincreasing. Some properties of monotone CAs have been discussed in [T2] . 

If a CA rule is number-nondecreasing and number-nonincreasing, it is 
number- conserving. The class of number-conserving rules is therefore a sub- 
class of the class of monotone rules. 

It is clear that all monotone rules are ENC, but have we a theorem stating 
a necessary and sufficient condition for a CA rule / to be monotone? The 
answer is "yes." It can be shown ^3], that 1 

A one-dimensional q-state n-input CA rule f is number-nondecreasing 
(resp. number-nonincreasing ) if, for all cyclic configurations of lengths 2n — 2 
and 2n — 1 it satisfies relation O (resp. 

Remark. The word and is essential since the CA rules that verify (j3J) or 
(jSJ) for all cyclic configurations of either length 2n — 2 or 2n — 1 are not 
necessarily monotone. To be monotone, say nondecreasing, a rule has to 
verify relation (J3J) for all configurations of lengths 2n — 2 and 2n — 1. This 
is, in particular, the case for three-state two-input CA rules. The sets of 
rules that satisfy relation (0J) for all cyclic configurations of lengths 2 and 3 

1 Anothcr, less practical, necessary and sufficient condition can be found in |12| . 



4 



have, respectively, 864 and 724 elements while the intersection of these two 
sets contains only 708 CA rules which are the only monotone nondecreasing 
rules. 

There exist 87 monotone elementary CA rules. The sets of nonincreasing 
and nondecreasing elementary CA rules both contain 46 rules. The intersec- 
tion of these two sets consists of the five number-conserving elementary CA 
rules. 2 

Although all monotone rules are ENC rules, are all ENC rules monotone? 
The answer is "no." In other words, the class of monotone rules is strictly 
included in the class of ENC rules. For example, elementary CA rules 99, 
173, and 229 are ENC but not monotone. They are defined, respectively, by 

/ 99 (0,0,0) = 1, / 99 (0,0,1) = 1, / 99 (0,1,0)=0, / 99 (0,1,1) = 0, 
/ 99 (1,0,0) = 0, / 99 (1,0,1) = 1, / 99 (1,1,0) = 1, / 99 (1,1,1) = 0, 



/ 173 (0,0,0) = 1, /i73(0,0,l 
/ 173 (1,0,0) = 0, / 173 (1 



0, / 173 (0,1,0) = 1, / 173 (0,1,1) = 1, 
• 173 (1,0,1) = 1, / 173 (1,1,0) = 0, / 173 (1,1,1) = 1, 



and 



/ 229 (0 , 0, 0) = 1, / 229 (0,0,1) = 0, / 229 (0, 1,0) = 1, / 229 (0,1,1) = 0, 
/ 229 (1,0, 0) = 0, / 229 (1,0,1) = 1, / 229 (1, 1,0) = 1, / 229 (1,1,1) = 1. 



3 Emulating number-conserving rules 

A CA rule f emulates rule g if, for all finite configuration of length L be- 
longing to the limit set A/ of rule f , its images by either f or g are equal. 
That is, for all c G Af and all i G Z L , 

g(c(i-r e ),c(i-r e + l), . . . ,c(i + r r )) = f(c(i-r e ),c(i-r e + l), . . .,c(i + r r )). 



It would be tempting to define a ENC CA rule as a CA rule that emulates, 
at least, one number-conserving rule. Here are two examples. 

2 See the appendix for a complete list. 
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Example 1. Elementary CA rule 176, defined by 

; 176 (o,o,o) = o, / 176 (o,o,i) = o, / 176 (o,i,o) = o, / 176 (o,i,i) = o, 

/ 17B (1,0,0) = 1, / 176 (1,0,1) = 1, / 17B (1,1,0) = 0, / 176 (1,1,1) = 1, 

is number-nonincreasing. One can readily verify that it emulates number- 
conserving elementary CA rules 184 and 240 (see top figure CJ. Configura- 
tions belonging to its limit set consist of isolated Is separated by sequences 
of 0s whose lengths depend upon the initial configuration. 
Example 2. Elementary CA rule 99 defined above, which is not monotone, 
emulates elementary CA rule 226. Configurations of its limit set consist of 
alternating sequences of 0s and Is with, depending upon the initial configu- 
ration, either a few pairs of 0s, separating two successive Is, or a few pairs 
of Is, separating two successive 0s (see bottom figure [TJ). 




Figure 1: The spatiotemporal patterns of elementary CA rules 176 (top) and 
99 (bottom) show clearly that these rules emulate, respectively, elementary 
CA rules 184 an d 240, and elementary rules 226. 

In Appendix 2 we show how to determine ENC CA rules emulating 
number-conserving rules. 

Are there ENC CA rules that do not emulate a number-conserving rule? 
Surprisingly, the answer is "y es ;" an d defining ENC CA rules as rules, which 
after a finite number of iterations emulate number-conserving rules would be 
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inadequate. As for monotone rules, the set of rules that emulate a number- 
conserving rules is strictly included in the set of ENC rules. For example, 
elementary CA rules 74 and 88, obtained by reflection 3 of rule 74, are exam- 
ples of such rules. They are defined by 

/ 74 (o,o,o) = o, ; 74 (o,o,i) = i, / 74 (o,i,o) = o, / 74 (o,i,i) = 1, 

/ 74 (1,0 , 0) = 0, / 74 (1,0,1) = 0, / 74 (1,1,0) = 1, / 74 (1,1,1) = 0, 

and 

/ 88 (0,0,0) = 0, / 88 (0,0,1) = 0, / 88 (0,1,0)= 0, / 88 (0,1,1) = 1, 

/ 88 (1,0 , 0) = 1, / 88 (1,0,1) = 0, / 88 (1,1,0) = 1, / 88 (1,1,1) = 0, 

They are not monotone and do not emulate number- conserving rules. Actu- 
ally, configurations belonging to their limit sets contains all the 8 different 
triplets. Both spatiotemporal patterns exhibit the propagation in opposite 
directions of similar structures (see figure EJ. Rules 173 and 229, defined 
above, and that are respectively conjugate of rules 74 and 88, do not emulate 
number-conserving rules and have spatiotemporal patterns that, like rules 74 
and 88, exhibit the propagation in opposite directions of similar structures. 

Since there exist ENC rules that are not monotone and do not emulate 
any number- conserving one, the union of the set of monotone rules and the 
set of rules emulating number- conserving rules is strictly included in the set 
of ENC rules. 



4 Motion representation 

The particle dynamics is not clearly exhibited by the CA rule table. A 
simpler and more visual picture of the evolution rule can be given by its 
motion representation. This representation, which has been first introduced 
for number- conserving rules |§J, was defined as follows. List all the neigh- 
borhoods of an occupied site represented by its site value s 6 Q. Then, for 
each neighborhood, indicate the displacements of the s particles by arrow(s) 

3 The g-state n-input CA rule R(f), denned by 

) = f(x n ,x n -i, ■ ■ -,Xi), 

is called the reflected of rule /. 
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Figure 2: Spatiotemporal patterns of elementary CA rules 74 (top) and 
88 (bottom). Both rules — one being the reflected rule of the other — are not 
monotone and do not emulate number- conservative rules. Initial configura- 
tions are random. 

joining the site where the particle (s) is (are) initially located to its (their) 
final position(s). A number above the arrow indicates how many particles 
are moving to the final position. To simplify the representation only neigh- 
borhoods for which, at least one particle is moving, are listed. Here are two 
examples. The motion representation of the two-state three-input rule 184, 
which represents the simplest car moving rule on a one-lane highway, is 

i 

rx 

10. 

Since a particle located on the left of an occupied site cannot move to the 
right (which the direction of motion) we do not mention it. In a less compact 
notation we could complete the above representation with 

o 

11. 

The motion representation of the three-state three-input rule 6171534259461 
is 

112 1 

rx rx rx rx 

10 11 20 21 
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In both examples particles move only to the right. It is, therefore, not nec- 
essary to indicate the state of the left neighboring site of the particle (s). 

It is rather straightforward to extend this representation to nonconserva- 
tive rules. The only difference with number-conserving rules is to indicate 
the possibility of creations and annihilations of particles. This is done by 
adding either a + or a — sign above the site state followed by a positive 
number n to represent, respectively, creation and annihilation of n particles 
at that site. Although the number of particles is not conserved, we shall keep 
the name motion representation, since creations and annihilations of parti- 
cles may be viewed as particles moving between the system and the external 
world. Here are two examples illustrating that ENC CAs rules can be viewed 
as evolving systems of particles and not just finite sequences of nonnegative 
integers obeying abstract transformation rules. 

Example 1. The motion of particles in a system evolving according to CA 
rule 176 is represented by 4 

A -i 
10 •, 01 

where the symbol • may represent either a or a 1. The second term shows 
that only isolated Is can survive, as illustrated in top figure [TJ In addition, 
when no more particles can be annihilated, that is, when all particles have 
empty neighboring sites, the evolution rule of the system coincides with CA 
rules 184 and 240, in agreement with the fact that rule 176 emulates these 
two rules. 

Example 2. The evolution of a system of particles evolving according to 
CA rule 74 is represented by 

001, (HO, 1 11. 

The first term shows that isolated Is move to the left (see top figure [2J, 
furthermore, starting from a triplet of Is in a sea of 0s, the combination 
of these three terms generates t a composite particle of length 6 oscillating 
between two states: 001110 and 011010 propagating to the left, or the two- 
row composite particle 

001110 
011010 

4 A systematic determination of motion representations is given in |15j . 



9 



which, then propagates to the left. Note that the motion representation lists 
only neighborhoods of moving particles, where "moving" stands for motion 
to or from the external world. 

5 Conclusion 

We have introduced a new class of CAs called eventually number- conserving 
(ENC) CAs adopting the following definition: A one- dimensional q-state n- 
input CA rule f is eventually number-conserving if, after a finite number 
of iterations, it satisfies the number- conserving condition for all cyclic 
configurations of length L > n. ENC CAs are models of open systems of 
interacting particles exchanging particles with the external world but that, 
after a finite transient time, the numbers of particles moving from and to 
the system exactly counterbalance each other. The interpretation in terms 
of particles is emphasized by the motion representation of the CA evolution 
rule. We have shown that the set of ENC CA rules strictly includes the 
set of number-conserving rules, the set of monotone rules, and the set of 
rules emulating a number-conserving rule. The set of ENC rules also strictly 
includes the union of all these sets showing that the set of ENC rules actually 
contains rules which do not belong to any of these sets. This challenging 
result suggests that the characterization of this new apparently simple class 
is probably not a trivial problem. 

Appendix 1 

The code numbers of monotone nonincreasing and nondecreasing elementary 
CA rules are respectively: 

{0, 2, 4, 8, 10, 12, 16, 24, 32, 34, 40, 42, 48, 56, 64, 66, 68, 72, 76, 80, 96, 98, 112, 
128, 130, 132, 136, 138, 140, 144, 152, 160, 162, 168, 170, 176, 184, 192, 194, 
196, 200, 204, 208, 224, 226, 240}, 

and 

{170, 171, 174, 175, 184, 185, 186, 187, 188, 189, 190, 191, 204, 205, 206, 207, 
220, 221, 222, 223, 226, 227, 230, 231, 234, 235, 236, 237, 238, 239, 240, 241, 
242, 243, 244, 245, 246, 247, 248, 249, 250, 251, 252, 253, 254, 255}. 
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There are five number- conserving elementary CA rules given by the intersec- 
tion of these two sets. Their code numbers are: {170, 184, 204, 226, 240}. 

Appendix 2 

In this appendix we show how to find examples of ENC CA rules emulating 
given number- conserving CA rules. 

Four-input rule 50358 is number-conserving. It is defined by 



/50358(0 , , , 0) 


= o, 


/ 5 0358(0 , , 0,1) 


= 1, 


750358 (0 , 0,1,0) 


= 1, 


/ 503 58(0, 0,1,1) 


= 0, 


/ 5 0358(0, 1,0 , 0) 


= 1, 


750358 ( 0,1,0,1) 


= 1, 


/ 5 0358(0,1,1,0) 


= 0, 


/ 5 0358(0, 1,1,1) 


= 1, 


750358(1,0 , , 0) 


= o, 


/ 5 0358(1,0,0,1) 


= 0, 


/ 5 0358(1,0,1,0) 


= 1, 


750358(1,0,1,1) 


= 0, 


/B035 8 (l, 1,0,0) 


= 0, 


/ 50 358(1,1,0,1) 


= o, 


750358(1,1,1,0) 


= 1, 






750358(1,1,1,1) 


= 1. 







As shown in figured its limit set, when the left and right radii are respectively 
equal to 3 and 0, can be viewed as 2-row tiles of one of the following types: 

1110000 1011000 
1011000 1110000 

concatenated with any number of 2-row tiles of type: 






Figure 3: Spatiotemporal pattern showing the limit set of J^-input rule 50358, 
where left and right radii have been chosen equal, respectively, to 3 and 0. 
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Among the 16 different quadruplets of Os and Is, the limit set of rule 50358 
contains only the following ones: 

(0,0,0,0), (0,0,0,1), (0,0,1,0), (0,0,1,1), (0,1,0,1), 
(0,1,1,0), (0,1,1,1), (1,0,0,0), (1,0,1,1), (1,1,0,0), 

(1,1,1,0). 

Among all the four-input rules having the same images as rule 50358 for all 
these 11 quadruplets one can verify that only rules 16566, 49334, and 50870 
are ENC and have the same limit set as rule 50358. These three rules are 
not monotone. 
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